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Abstract 

A new formalism using the ideas of the augmented space recursion (introduced by 
one of us) has been proposed to study the ground state properties of ordered and dis- 
ordered one-dimensional Holstein model. For ordered case our method works equally 
well in all parametric regime and matches with the existing exact diagonalization 
and DMRG results. On the other hand the quenched substitutionally disordered 
model works in low and intermediate regime of electron phonon coupling. Effect of 
phononic and substitutional disorder are treated on equal footing. 

Key words: Holstein model, Polarons 
1 Introduction 

When a conduction electron or hole moves through a polar crystal it polarizes and 
distorts the neighboring ion lattice sites. When it moves it carries along with it the lattice 
distortion. The electron together with this lattice distortion or self-consistent polarization 
field is a quasi particle and is referred to as polaron. Many experiments indicate strong 
evidence of polaronic carriers in strongly correlated electronic materials like manganites 
(showing colossal magneto-resistance) [l]-[2], organic [3] and high Tc superconducting 
cuprates [4] -[5]. Holstein model [6] has been one of the basic model used to study the 
formation and nature of polaron. Analytical approaches to solve the model have got 
limited success. Main draw back of those results are they are valid only in a limited range 
of parametric regime of electron-phonon coupling [7] and also for small system size. 

Recently there has been a spurt in numerical studies on this model, which includes 
the variational approach based on exact diagonalization method (VAED) [8]- [9], the 
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density-matrix rcnormalization group techniques (DMRG) [12], exact diagonalization 
techniques(ED) [13]-[19], the quantum Monte Carlo calculations (QMC) [20]-[21] and 
global- local method [22]. None of these methods work equally well in all parametric 
regimes. Each has some shortcomings. VAED, the most recent and successful of the meth- 
ods, though highly accurate in the the weak and intermediate coupling regimes, fails to 
maintain its accuracy in the strong coupling limit. The DMRG method is accurate in the 
strong and intermediate coupling regimes, but is not accurate in the weak coupling regime. 
The DMRG results arc not translationally symmetric. DMRG lacks the accuracy of VAED 
in the weak coupling regime because its lattice size is smaller than the spatial extent of 
the large polaron that will be formed. The VAED cannot emulate its success of the low 
and intermediate regime in the the strong coupling regime because the electron-phonon 
system demands more phonons in its 'root state' [9] and its near vicinity. 

Disordered Holstein model has been earlier attempted by [10] using Dynamical CPA. But 
their method is effective in the intermediate electron-phonon coupling regime and for only 
infinite coordination, i.e. on a Bethe Lattice. Keeping this in mind, in this communication 
we have proposed a reciprocal space recursion technique, introduced earlier by us in a 
different context [11], as a feasible method for the study of the Holstein model. We have 
constructed a basis which meets the requirement in all parametric regimes. 



2 The Augmented Space Formalism for ordered Holstein model 

In this section we shall describe the setting up of the electron-phonon Hamiltonian as 
an operator on the underlying Hilbert space. The electronic part will be represented in 
a tight-binding like basis. In our simple model there will be only a single orbital per 
site and the tight-binding basis will be labeled by a site i?„ : This basis will 

span the electronic part of the Hilbert space H^'. The phonon part will be described by a 
configuration state which will indicate how many phonon excitations there are at each site 
labeled by Rn- Each conGguration will be a pattern of the type {tir^ , ur^ . . . ur^ . . .} where 
ur^ takes the values 0, 1, 2, . . . and is called the cardinality of the site The pattern 
is called the cardinality sequence and uniquely labels a phonon state. These states will 
span the phonon part of the Hilbert space H^^. The total Hilbert space of states of the 
electron-phonon system will be H = H^' (8) H^'*. A member of the basis in this product 
space will be labeled both by the site labeled real-space part and a configuration part 
which tells us how many local phononic excitations are there at each site in the lattice, 
e.g. \Rn ® W,,nR2, . . .}). 

In this basis, the Holstein Hamiltonian [6] is : 

eY,Pn®I-tY.Y.Tnm® I + hiJ^ I ^ - Pn ^ (T+ + (1) 

Rri Rm Rn Rn 

where, 

• Pn is the projection operator |i?„)(i?„| and T„„ = \Rn){Rm\ is the transfer operator in 
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the electronic part of the Hilbert space H*^', 

• Tj^ = I • • • i^Rn ± 1 •••)(•• • nRr, • • • I are the step up or step down operators and Nn is the 
number operator riji^ \ . . . tir^ ...){... tir^ . . . | both in the phononic part of the Hilbert 
space H^'*. 

• Here a single electron with an on-site energy e can hop to its near neighbour on an in- 
finite lattice with hopping integral t and interacts with dispersion less optical phonons 
with frequency u. The clectron-phonon coupling strength is denoted by A . The param- 
eters e,t, hw and A have the units of energy 

• In all subsequent calculations energies are scaled in units of t. 

This type of Hamiltonian in an augmented space H^' ® "HP^ has been introduced earlier 
by us [23] to deal with time dependent disorder in an electronic system, precisely of the 
form produced by lattice vibrations. The great advantage of this formulation is that if we 
have quenched disorder in e^Rn)-, then a simple augmentation of the Hilbert space H^' by 
the configuration space of the e{Rn) allows us to deal with all the configuration averaged 
properties of a disordered Holstein model. This extension will be carried out by us in a 
later section. Since we are interested in calculating the spectral function, dispersion curves 
of the system, we generalize the above mentioned formalism to generate basis states in 
reciprocal space. Here a general state in the augmented reciprocal-space basis has the 
following form : 

\k ® {uR^.nR^ ...ur^...})^ (l/y]V) ^exp(iA; • Rn)\Rn ® {ur^^ur^ ...ur^...}). 

Rn 

Here /c is a reciprocal space vector and we have Fourier transformed the electronic part 
of the basis to momentum space keeping the bosonic part unchanged. In the recursion 
method, which we shall later employ, we shall generate the basis states by repeatedly 
applying the Hamiltonian on the following starting state : 

|1) = \k® {0}) = (l/ViV) ^exp{^A; • Rn]\Rn ® {0}) 

Rn 

Here {0} denotes a zero phonon configuration in the bosonic subspace of the full Hamil- 
tonian. That is, with a null cardinality sequence. 

Alternatively, we can generate basis states recursively by similar action of the Hamilto- 
nian on the following two states : 

(i) |1) = \k® {0}) = {l/^/N)Y,^^v{^k ■ Rn}\Rn ® {0}) 

Rn 

(ii) |2) = |A; (8) {TVi, 0, ...» = (1/v^) ^ exp{ik ■ (g) {N,, 0, . . .}) 

Rn 

Here A'^1,0,0.. denotes A^i number of phonon at site 1 of the bosonic subspace of the 
Hamiltonian and zero else where. We shall call this the phonon-enrichcd state. 

The advantage of the second choice is that, if we apply the Hamiltonian N times on the 
two initial states, -|- A^i number of phonons are generated in the site 1 and sufficient 
number of bosons are generated in its vicinity, while in the first choice we get only N 
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number of phonons in the first site and progressively less numbers in the neighboring 
sites. This is important in order to get a convergent result in the high e-ph coupling 
regime. 

In order to calculate averaged spectral function we tri-diagonalize the Hamiltonian using 
the following three term recursion : 

We start from a state \k ® {ur-^, . . .}) = 1 1) and generate the other states through the 
following three term recursion : 

|n + 1) = H\n) - an(k)\n) - Pl_,(k)\n - 1) (2) 

The coefficients an{k) and P^i^) obtained by ensuring that the new basis generated 
is orthogonal. This yields : 

<^n{k) = {n\H\n) / {n\n) Pni^) ~ + + 1)/ {n\n) 

To carry out the recursion we need to know the action of the Hamiltonian (1) on a 
general state : If write the Hamiltonian (1) as : H — Hi + H2 + H3 + H4 then 

Hi\k (g) {ur^^ur^ . . .])^e\k ® {uR^.nR,^ . . .}) 

H2\k ® {ur^^ur^ . . .}) = -t^exp{iA; ■ x}\k 'S) {ur^^-^^ur^^^ . . .}) 

X 

H3\k® {nR„nR^...})= ^^X)?^Rnj \k ® {nR^,nR^ . . .}) 

H4\k {nR^,nR^ . . .}) = -A|A; (g) {ur^ ± 1,71^2, . . .}) (3) 

Once we carry out a recursive determination of the coefficients {an{k) , (3^{k)} , the Green 
function is given as a continued fraction : 

G{k,E) = ' 



m 



E-a,ik) m. 



E - aj,{k) - /3i,{k) T{k,E) 

We carry out recursion up to a finite number of N steps and then terminate the continued 
fraction with a herglotz terminator T{k, E) [26] suggested by Beer and Pettifor [27]. The 
spectral function is then given by : 

A{k,E) = -(l/7r)!3m < G{k,E) > 
The dispersion curves are obtained by fitting Lorentzian's in the neighbourhood of the 
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peaks in the spectral function : the Lorentzian centre gives the energy E(k) while the 
width gives the width in the dispersion curves. 

We apply the recursion based conjugate gradient (CG) technique [25] to obtain ground 
state energy and wave function. To calculate the excited state wave function we find a 
orthogonal state to the ground state using Gram-Schmidt method of orthogonalization 
and then apply CG technique to find the first excited state energy and wave function. 
With the wave function and dispersion curve at our disposal we find the various correlation 
functions, effective mass from it. The effective mass we have calculated using the standard 
formula : 



mo 1 d'^E{k) 



m* 2t dk^ 



(4) 

fe=0 



The mean phonon number indicates a measure of the phononic character of the polaron. 
N^'^ =(V'o|/«)A^n|V'o> (5) 

Here |^o) is the ground state wave function obtained from the CG technique. The static 
correlation function between the electron position and oscillator displacement is given 

x{Ri - Rj) = {HP^ ® {T; + T-)\iJo) (6) 

The number of excited phonons in the vicinity of the electron is given by 

^(Ri-Rj) = {iljo\Pi®Nj\iljo) (7) 



3 Augmented space formalism for disordered Holstein model 

In this section we shall generalize our ideas to a disordered Holstein model. We shall 
start with a Hamiltonian : 

Rn Rn Rm Rn Rn 

where, 

£{Rn) = ea n{Rn) +83 {1- n{Rn)) 

Here n{Rn) is a random variable which takes the value if Rn is occupied by an A type 
of atom and 1 if it is occupied by a B type. The probability of these events are x (the 
concentration of A atoms in the alloy) and 1 — x (concentration of B atoms in the alloy) 
respectively. IN this simple model we shall consider only this type of binary substitutional 
disorder. For real disordered alloys the other parameters t, cu and A may also be random. 
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The augmented space formalism for quenched disorder [28] writes the probabihty density 
of the random variables as : 



P[n{Rn)] =x 6[n{Rn) - 1] + {1 - x) 6[n{Rn)] 

--(I/tt) S^m (T I ({n{R^)+zS)I-M(^A''\]) 



The operator M^") associated with the random variable n{Rji) is such that its spectral 
density is the probability density of the random variable. Here the representation of this 
operator is 




in the basis | j) = v^|0) + Vl ~ and | J,) = ^/l — x\0) — \/x\l), where |0) and |1) are 
the eigenvectors of M^") with eigenvalues and 1 respectively. 

The configuration space of a single random variable n{Rn) : ^"-"^ is of rank 2 and is 
spanned by these vectors | t) and | |) and M^") is an operator on this space. We consider 
the configuration space of all the variables {n{Rn)} : $ = n® 0^"^- The augmented space 
formalism [28] constructs a Hamiltonian : 



Rn Fin Ftn 

-^EE^nm® /®/ + ;za;^/®iV„®/- A^P„® (t+ + T„) ®/ (8) 

Rn Rm Rn Rn 

where, 

• -Pf = I I and T|" = | t")(i"' | + | i"')(T"' | are the projection and transfer operators 
in the configuration space $ 

• ei = (1 - x){eA - £b) and 62 = \Jx{l - x){eA - £b)- 

This enlarged Hamiltonian is in the full augmented Hilbert space : H = H"^' (E) ® 
As in the case of the phonon space, in this disorder configuration space, a general vector 
is a pattern of t and J,-s. The sequence of sites where we have a J, is called the cardinahty 
sequence and uniquely describes the configuration. A typical member of a basis in this 
augmented space is \Rj (8) {n^^, n^j, . . .} (8) {Rki,Rk2J ■■■})■ 

The augmented space theorem [28] states that : 

< G{R, E) >= {R ® {0} ® 0| (ei -Hy^ |i? ® {0} ® 0) (9) 

where |0) e $ is the null cardinality sequence or one where we have t everywhere. Equiva- 
lently, exactly as discussed in section 1, we can construct the above formalism in reciprocal 
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space by Fourier transforming the electronic part of the basis and keeping the rest un- 
changed. So we can write the configuration averaged Green's function in reciprocal space 
as : 

< ^(A;, £;) >= (A; {0} 0| {eI -H)'^ |A;(8){0}®0) (10) 



The operation of the terms in the Hamiltonian are mostly identical to (7) except for the 
second to fourth terms in equation (8) : 



H2\k ® . . .} ® {C}) = £i 5{Ri e {C}) \k ® {UR,,. . .} ^ {C}) 

H3\k^{nR^,...}(^{C}) = 82 \k®{nR^,...}®{C±R^}) 

Hi\k ® {hr^, . . .] ® {C}) = -tY,Q-^v{ik-x]\k®{nR,-X---)®{C-x]) 



(11) 



here, {C} is a general cardinality sequence {Rk, Rm, ■ ■ ■}■ 

For obtaining the configuration averaged spectral function we carry out recursion with 
the starting state \k {ur^, . . .} (8) 0) = |1) using the three term recursion similar to 
(6) but in the full augmented space. As before, from the orthogonality of the recursive 
states, we obtain the coefficients a{k),P{k) and averaged Green function is obtained as a 
continued fraction : 

« ^(^' ^) »= (12) 



E - a2{k) - 



E-aN-(3i,T{k,E) 



As before the calculation is carried down to a finite number N of steps and then the 
continued fraction is terminated by a function T(A;, E) as suggested by Beer and Pcttifor. 
The configuration averaged spectral function and other quantities of interest are calculated 
using the same formalism as discussed in section 1. 



4 The Ground state in a ordered Holstein model 



In this section we will discuss the results of our ground-state properties and compare 
them with the results of two of the most successful recent numerical works. 

The figure 1 shows the spectral functions for different k values. All the spectral functions 
show a very narrow delta function like peaks at the lower end of the spectrum. This is 
related to the band corresponding to the polaron ground state. This state has very narrow 
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Energy(E) 



Fig. 1. The spectral functions are plotted for different k values. The spectral functions are 
obtained from a k-space recursion with a 22 shell nearest neighbour map. 

width, which means this fc-labeled state has a very long lifetime. The other peaks are 
related to excited states. These have larger widths. 

The dispersion curves shown in figure 2 (left) are obtained by fitting Lorentzian's in 
the neighbourhood of the peaks in the spectral function : the Lorentzian centre gives the 
energy E{k) while the width gives the width in the dispersion curves. 

Table-I compares our polaron ground-state energy at k=0 with those obtained by VAED 
[8] (in the weak couphng regime) and DMRG [12] (in the Strong coupling regime) for two 
sets of parameters. The ground state energies were obtained by two different starting 
root states as mentioned earlier. Our energies match extremely well with both the earlier 
VAED and DMRG calculations. 

We have calculated the effective mass spanning all parametric regimes. We find good 
agreement with results of the DMRG and the earlier VAED methods. This is shown in 
figure 3 (left panel). To check numerical stability we have done the calculations first with 
a 17 shell map with a maximum of 34 phonons and another with a 18 shell map with a 
maximum of 35 phonons. The results match to within our accuracy window across the 
parametric regime < g < 5. At the low electron-phonon coupling regime we have a 
quasi-free electron with a slightly renormalized mass. As the electron-phonon coupling 
strength increases the polaron becomes heavier. The crossing though always smooth, is 



X/u 


Present (VAED) 


Present (M- VAED) 


VAED [8] 


DMRG [12] 


1.0 


-2.40968172.3933 


-2. 1G9G8 1723933 


-2. 1G9G8 1723933 


-2.1G9G8 


V2 


-2.998828186866 


-2.998828186866 


-2.998828186867 


-2.99883 



Table 1 

Comparison of Polaron Ground state energy for k=0 for two different parameters. M-VAED 
refers to the phonon enriched starting state described in section 1 
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1=0.5 



E(k) 



■2.5 
■2.7 
■3 
■3.3 

■9.11 
-9.12 

■9.13, 



g=1.0 



g=l.5 



g=3 



0.2 



0.4 



0.6 



0.8 



Fig. 2. (left) Dispersion curves for ground state and few excited states are plotted for A = a; = 1. 
The recursion was carried out using a 22 shell map. (right) Dispersion curves for different values 
of e-ph coupling. Here X/u = g 




Fig. 3. The left panel shows the logarithm of the ratio of the polaronic mass m* to bare electron 
mass w.r.t. the electron-phonon coupling strength in different Limits. Dots are the DMRG 
extracted data points and Squares are VAED extracted data points. The right panel shows the 
average kinetic energy for different values of lu 

rapid in the adiabatic limit. In the same figure (right panel) we also plot the averaged 
kinetic energy, which rapidly decreases as the polaron becomes heavier . 

We have calculated the static correlation function x{Ri ~ Rj) between the electron po- 
sition and quanta of lattice vibration which gives us a measure of the electron-induced 
lattice deformation and its spatial extent. We compare our results with VAED and DMRG 
results. There is excellent agreement in between our results generated using both of our 
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basis and VAED results which is shown in figure 4. In figure 4(a) we show the static 
correlation function for uj=0.1 (adiabatic) and A=0.1(weak coupling), and we can clearly 
see spatial extent of the polaron span the system size, therefore a large polaron. We do 
not show the DMRG result in this case as it is less rehable in this parameter regime. In 
figure4(b) we show the x for co—l.O and A=0.5, which again matches well with the VAED 
calculation. Figure4(c) and figure4(d) give us the result for high electron-phonon coupling 
in adiabatic uj/t=0.1 (A=4.35) and intermediate cj/t=1.0(A=3.0) limits respectively. In 
both the cases spatial extent of the polaron has been reduced, resulting in small polaron. 
Here we have achieved excellent convergence applying the modified or phonon-enriched 
VAED basis (M-VAED), as described in section 1, using N=17 (34 bosons at the root 
site); N=18 (36 phonons at root site) and N=19 (37 phonons at the root site) shell maps. 
Our energy converges to 11-12 decimal places for these two calculations. Figure4(c) shows 
good agreement with DMRG results. The value of our extracted DMRG x(0) is 5.459016. 
The value of our x(0) for this case is 5.374, which is very close to 5.4 the lower hmit of 
extrapolated VAED data. Figure4(d) also has a good agreement with DMRG results. The 
value of our x(0) is 5.91 and value of the DMRG extracted data is 5.73. For small polaron 
our data are in good agreement with DMRG data. Unlike DMRG results , our results are 
refiection symmetric (i.e.,x(l)=x(-l)) as this basis takes proper care of translation sym- 
metry. We conclude that we have achieved highly reliable results in all regimes. 

We have also calculated x at non-zero k at a;=0.8, A^=0.4 (figure 5, left) and a;=1.0, 
A=3.0(figure 5,right) to see how the polaron transforms from predominantly electronic 
character at k=0 to phononic at k=7r. In figure 5(left) at k=0 where group velocity is zero 
the deformation affects only the sites in its close vicinity , falls off exponentially and is 
always positive. At k=7r/4, there is a enhancement in the deformation amplitude and it 
acquires a negative sign in oscillation as well. At k=7r/2 the oscillations are enhanced in 
sign as well as spatial extent. The spatial extent of the deforming oscillations attains its 
maximum at k—ir. Our calculation matches well with VAED for same set of parameters 




0.4 
0.3 
0.2 
0.1 




-10 





(b) 


M-VAED , 




VAED-P / 




VAED / 













Nh=17 




Nh=18 1 




DMRG 1 




■ J 





'10 



(i-i) 



-10 



10 



--- Nh=17 




Nh=18 j 




o DMRG 


1 " 







10 



Fig. 4. The Lattice deformation x as function of Ri 
and A=.5, (c)a;=.l andA=.435, (d)a;=l. and A=3. 



Rj at k=0 for (a) uj=.l andA=.l, (b)ti;=l. 
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(H) (i-i) (i-i) (i-i) 

Fig. 5. (left) Figure shows lattice deformation x foi' different k values at cj = 0.8 and = 0.4. 
(right) Figure shows lattice deformation x for different k values at w = 1.0 and A = 3.0 



for zero as well as non-zero k-values. 

Figure 5 (right) shows the x same k-values but for high electron-phonon coupling 
strength. Here x shows hardly any variation with k, i.e., lattice distortion is predominantly 
on the electron site, with hardly any group velocity. Figure 2 (right panel) shows polaron 
is almost dispersion less for this parameter i.e., group velocity((ia;/(iA;) is almost zero and 
again the average kinetic energy(figure 3 right panel) also drops drastically for this value, 
hinting clearly that a very small polaron has formed. 

We have also calculated the correlation functions x(0), xiX)-, (x(0)-x(l))/g which is 
shown in figure 6. Where g is A/cu. x(0) has a non-linear behavior in the intermediate 
electron-phonon coupling regime. It is linear w.r.t g both in the low and high electron- 
phonon coupling regime and this trend is there in all limits(i.e., for different uo). But slope 
of x(0) in the low electron-phonon coupling regime depends on a; but its slope in high 
electron-phonon coupling regime becomes independent of uj and its value approaches 2. 
The change in slope of x(0) as we go from low to high g, is more prominent for lower 
values of UJ. x(l) increases linearly to start with(i.e., low g regime), then its rate of 
increment decreases and then it starts decreasing, signaling the arrival of high electron- 
phonon coupling regime. The decrement in %(!) w.r.t to g and the corresponding change 
in x(0) imphes that, from this g onward the lattice distortion starts getting confined to the 
electron site. The correlation function (x(0)-x(l))/g substantiates these facts. It always 
saturates to a value 2 for high g regime. For c<j=0.1 we have not gone to that magnitude of 
g where (x(0)-x(l))/g would saturate to 2, but the trend is quite clear. For lower values of 
UJ the cross-over is rather sharp, which is also reflected in the Effective mass and average 
kinetic energy calculation (figure 3). 
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co=0.1 a>=1.0 co=4.0 




X/(o X/d) X/o) 



Fig. 6. Figures in the upper row shows the x(0) at three different 00. Middle row is for x(l). 
Lower row is for (x(0)-x(l))/g 

5 The First Excited-State in an ordered Holstein model 

In this section we discuss our calculation for the first excited state of the electron-phonon 
system. The first excited state consists of the ground state polaron and an unbounded 
extra phonon excitation [8]. Here the energy difference between the ground state and the 
first excited state should be equal to uj for a infinite chain(i.e., thermodynamic limit) 
and the mean phonon number difference should be equal to one(i.e., AN^^^''=1.0). But for 
the excited state there are two distinct regimes, one below the critical electron-phonon 
coupling Ac where the excited phonon tries to be at a infinite distance away from the 
ground state polaron(but the finite size of the system hinders it)and above it where 
the excited phonon is absorbed by the ground state polaron forming a bound state. The 
calculations here are done for a;=0.5. VAED[8] calculations showed that a phase-transition 
occurs in the first excited state at Ac=0.95(i.e., g=1.9) for a;=0.5. Our calculations shows a 
phase-transition precisely at A=0.95. Our calculated binding energy ( for the first excited 
state) and correlation functions clearly answers to the issue, whether the extra phonon 
excitation forms a bound state with the ground state polaron after a certain point (i,e., Ac) 
or prefers to remain infinitely separated (of course limited by the finite size of the system) 
for all values of A. 

The binding energy A = AE - (where AE = Ei — Eq, Ei and Eq are the first-excited- 
state and ground-state energy respectively) as a function of electron-phonon coupling A 
is shown in figure? for various system sizes. Below Ac, A varies with system size but is 
greater than zero. This variation with system size is due to the fact, that the excited 
phonon want to be infinitely away from the the ground state polaron, but the system size 
limits its ability. As the system size increases it slowly approaches the thermodynamic 
limit, i.e., A=0. For A > Ac, where the absorption of the excited phonon by the ground- 
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K 

Fig. 7. Figures shows the binding energy as a function of A for different system sizes(i.e., using 
different shell maps). 




X/at=g X/(o=g 

Fig. 8. Figures show the comparison between the first excited state and ground state Correlation 
Functions at a;=0.5. Phase transition in the first excited state occurs at 7=0.95 i.e., g=1.9. Solid 
line with dots are results for the first excited state and dashed line with squares are results for 
the ground state. 

state-polaron has resulted in formation of a bound state, A has clearly converged at 
Nh=14 and is negative. 

Figures compares the different correlation functions of the first excited state and the 
ground state for a;=0.5. There is steady difference in between N^^ of the first excited and 
the ground sta,te{ANP^ ~1.0 ) below g=1.9, and then starts diverging. They again appears 
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Fig. 9. Figures shows the phonon number 7 as a function of the distance from the electron 
position i?j — Rj . 

to converge at higher g though our excited state are not very accurate for higher g. Average 
kinetic energy, x(0), x(l) are same for the ground state and the excited state below Ac 
which very clearly indicates that below Ac the excited phonon and the ground state polaron 
are unbounded. The average kinetic energy, x{^): of the first excited state below Ac 
are entirely that of the ground state polaron and it is transparent to the presence of the 
excited phonon[8]. At Ac the root site disassociates itself from the rest of the lattice[8], 
the signature of which can be found in the sudden rise in x(0) and corresponding fall in 
. The bound state formed by the absorption of the excited phonon by the ground 
state polaron is a excited polaron. It slowly stabilizes with increasing A and exhibits the 
behavior of a de-excited polaron. 

In figure 9 we have calculated the distribution of the number of excited phonons in the 
vicinity of the electron '-/{Ri — Rj). We show '~f{Ri — Rj) both for ground state (70) and and 
excited state(7i) just above(A=1.0) and below (A=0.9) the transition point. Below the 
transition point the peaks of 70 and 71 are almost same , but 71 has a longer tail suggesting 
that tail represent the extra excited phonon which extends throughout the system in its 
attempt to remain unbounded from the ground state polaron. Our calculated Nph are 
same as VAED[8] for all the cases. The difference in between the ground state Nph and 
first excited state Nph should be one for a infinite system, but our difference is about 1.02 
and this can be attributed to finite size effect [8]. The situation above Ac is different. Here 
the peak value of 71 is almost three times the peak value of 70 and secondly here 71 decays 
fast unlike that below Ac- Here the difference in phonon number of the ground state and 
the excited state is 2.33 again exactly same as VAED[8]. Thus here the excited bound 
polaron has many extra phonon excitations compared to the ground state polaron. Since 
the extra phonon excitation are almost confined to the root site of the electron, it throws 
some light on the fact that at Ac the root site gets detached from the rest of the lattice[8] . 
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6 Results for a disordered Holstein model 




Fig. 10. The spectral functions for the disordered Holstein model 
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Fig. 11. Comparison of the average phonon number and kinetic energy for the disordered and 
ordered Holstein model at different oscillator frequencies 
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Fig. 12. (Top) The dispersion curves for the disordered Holstein model, for g=l, e^=0, 6^=0.5, 
for three different compositions. (Bottom) The dispersion curves for the disordered Holstein 
model, for g=0.01, €^=0, 6^=0.5, for the 50-50 alloy. 

In figure 10 we show the spectral function from the F to X point in Brillouin zone 
for the 50-50 alloy. The spectral functions show the extra disorder induced widths as 
well as contributions coming from both the components. This is much better seen in the 
dispersion curves shown in figure 12 (top). The dispersion curves for both the disordered 
alloys show the two branches arising out of the two components and the large disorder 
induced widths in the upper band. In the low concentration impurity regime the two 
bands almost cross each other. For higher concentrations the bands are well separated. 
As in the ordered case, the lower band has very little width. In comparison with the 
ordered case the upper band has large width particularly near the Brillouin zone edge. 
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co=0.1 co=1.0 co=4.0 




Fig. 13. Comparison of x(0) and x(l) for the disordered and ordered Holstein model at different 
oscillator frequencies 

This extra width arises from quenched disorder scattering. These dispersion curves are 
for g=l. Figure 12 (bottom) shows the 

dispersion curves for g=0.01, i.e. low electron-phonon coupling. Here the disorder ef- 
fect dominates and both the branches show disorder induced large smearing. For low 
electron-phonon coupling in a disordered alloy, the polaron, when formed, will have dis- 
order scattering induced finite life-time. We conclude that disorder scattering effects are 
small on heavy polarons, while they tend to give larger lifetime effects to light polarons. 

We have also calculated different ground state correlations functions for the disordered 
Holstein model and have compared it with its ordered counterparts to get a better insight 
of the disorder effect. Figure 11 compares the average phonon number and the average ki- 
netic energy for the ordered and disordered case for three oscillator frequencies. Site disor- 
der tends enhances the average phonon number with increasing electron-phonon coupling 
in the adiabatic regime where as, in the intermediate and anti-adiabatic regime disorder 
has hardly any effect. Disorder lowers the average kinetic energy in all the regimes. Figure 
13 compares the x(0) and x(l) and similar trend is observed in this' set of correlation 
functions too. We conclude that the disordered effect is more prominent in the adia- 
batic regime and with increase in the oscillator frequency, the phononic disorder becomes 
dominant. 
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